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ABSTRACT 

The trajectories of Solar Energetic Particles (SEPs) in an Interplanetary Magnetic Field 
(IMF) exhibiting large-scale fluctuations due to footpoint motions originating in the photosphere, 
are simulated using a full-orbit test-particle code. The cross-field transport experienced by the 
particles in three propagation conditions (scatter-free, with scattering mean free path A=0.3 AU 
and A=2 AU) is characterized in the Parker spiral geometry. The role of expansion of the magnetic 
field with radial distance from the Sun is taken into consideration in the calculation of particle 
displacements and diffusion coefficients from the output of the simulations. It is found that 
transport across the magnetic field is enhanced in the A=0.3 AU and A=2 AU cases, compared to 
the scatter-free case. Values of the ratios of perpendicular to parallel diffusion coefficients vary 
between 0.01 and 0.08. The ratio of latitudinal to longitudinal diffusion coefficient perpendicular 
to the magnetic field is typically 0.2, suggesting that transport in latitude may be less efficient. 



Sun:particle emission, Sun:heliosphere, Solar-terrestrial relations, Solar wind, 



Subject headings: 
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Introduction 



As Solar Energetic Particles are transported 
through the Interplanetary Magnetic Field out 
to locations at 1 AU and beyond, they experi- 
ence magnetic field fluctuations which occur over a 
range of scales. The main consequences of a parti- 
cle's interaction with such fluctuations are thought 
to be: (a) diffusion in the direction parallel to the 
magnetic field, characterized by a diffusion coeffi- 
cient km, that typically describes pitch-angle scat- 
tering counteracting the strong magnetic focusing 
in the IMF and (b) diffusion in the direction per- 
pendicular to the average field, characterized by a 
coefficient n± , thought to be orders of magnitude 
smaller than K|| . 

Diffusion in the parallel direction has been the 
subject of a large number of studies over the past 
three decades, relevant not only to SEPs but also 
cosmic rays and more general particle propaga- 
tion in turbulent fields. While the initial consensus 



was that diffusion was the main factor determining 
profiles of SEP intensity versus time, a paradigm 
shift took place, that implied the role of interplan- 
etary traiisporHn shaping SEP profiles was mini- 
mal (|Reameslll99"9h . More recent studies, however, 
have re-affirmed that diff usion does play a ro le in 



determining SEP profiles (jMason et al. 2006) 



For SEPs, diffusion in the perpendicular direc- 
tion has been largely ignored until recently, and is 
currently the subject of much interest and contro- 
versy. Some of the observational studies pointing 
towards perpendicular transport playing an im- 
portant role, used data at h i gh heliolatitudes fr om 



Ulysses (| Zhang et al 
Recent observations 
electrons in 



2009t balla et all l2003allbh 



from STEREO show that 
3 He-rich events reach spacecraft 
widely separated in longitude, a fact that may re- 
quire transport across the field (jWiedenbeck et aT 
20101) . A number of theoretical studies have also 



addressed perpendicular diffusion for SEPs, as de- 
tailed below. 
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A common technique when modelling energetic 
particle propagation is to adopt a kinetic approach 
with a diffusion tensor including parallel and per- 
pendicular diffusion coefficients, Ku and k±. In 
this approach Kit and k± are an input to the model 
and are estimated using observational data. Par- 
ticle gyromotion is neglected and a kinetic trans- 
port equation is solved numerically to give a gyro- 
averaged phase-space density function for the SEP 
population. 

Such an approach is taken by IZhang et al 



(j2009h where the Fokker-Planck equation is re- 
cast into a set of stochastic differential equations 
(SDEs) that are solved numerically to investigate 
the propagation of high-energy SEPs. They found 
that the inclusion of perpendicular diffusion in- 
creased the uniformity of particle flux when calcu- 
lated at different locations during the decay phase 
of a simulated event. This phenomenon has been 
observed in SEP events and is refered to as the 



"reservoir phenomenon" (|McKibbenlll972l) 



A similar approach is taken in iHe et al.l (|201lh 
where the influence of SEP source characteristics 
on observations at 1AU is investigated. They 
found that perpendicular diffusion plays a signifi- 
cant role in SEP propagation, particularly in cases 
where a spacecraft is not directly connected to the 
acceleration region; in such cases, the earliest ar- 
riving particles can be seen propagating towards 
the Sun, having scattered backwards at large dis- 
tances. 

(|2010tl 



Drbge etal 



solved the Fokker-Planck 
equation using a time-forward SDE method with 
two different particle scattering schemes, to inves- 
tigate the spatial distribution of SEPs in the inner 
heliosphere. They found that scaling the value of 
Aj^ - the mean free path in the direction perpen- 
dicular to the field - with gyroradius gave more 
realistic spatial distributions at 1 AU than when 
Aj_/A|| was kept constant at a value of 0.01. 

Another technique that can be used to under- 
stand the cross-field transport of SEPs is a full- 
orbit test-particle simulation. In this scheme a 
particle's equations of motion are solved numer- 
ically to determine its trajectory. Such an ap- 
proach assumes that test-particles do not interact 
with one another, or affect the IMF. 



port in a model of the IMF that included isotropic 
small-scale turbulence. They found that such 
turbulence led to diffusive behaviour and that 
Aj_/A|| = 0.1 at around 10AU, slightly larger than 
observed values. 

In this paper, the propagation of SEPs in large- 
scale turbulence is investigated by means of full- 
orbit test-particle simulations. The turbulence is 
induced by footpoint motion associated with su- 
pergranulation, leading to large-scale fluctuations, 
according to the description first developed by 
iGiacalond (|1999l ) to investigate the propagation 
of CIR-associated particles to high heliolatitudes. 
He ran full-orbit test-particle simulations of 1 McV 
protons, injected at 2 AU and integrated for a pe- 
riod of 10 days, and found that field line braiding 
resulting from magnetic footpoint motion provides 
an efficient mechanism for latitudinal transport. 



Pei et al.l (|2006f ) used a test-particle method to 



investigate the onset ti mes of SEPs in the IMF 
model put forward by iGiacalone ( 19991 ). They 
found that such fluctuations can reduce the length 
of magnetic field lines and hence lead to onset 
times for SEP events shorter than within a Parker 
spiral model. 

Here the effects of large-scale turbulence on 
the spatial distribution of SEPs, in particular 
its role spreading particles along the direction 
perpendicular to the magnetic field, are investi- 
gated. Three propagation regimes are analysed: a 
scatter-free regime and two scattering ones with 
different mean free paths. Diffusion coefficients 
are calculated from the output of the test parti- 
cle model and the influence of the expansion of 
Parker spiral magnetic field lines on their values 
and evolution discussed. 

An overview of the IMF model and the param- 
eters used in simulations is presented in Section 
[21 along with details relating to the test-particle 
method's implementation. In Section[3]the results 
of scatter-free simulations are presented and con- 
trasted with corresponding simulations which in- 
clude scattering. Finally in Section 0] the results 
are summarized and conclusions are outlined. 



Tautz et alj ([20111 ) used a full-orbit test-particle 



method to investigate the properties of SEP trans- 
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2. Model 



2.1. Interplanetary Magnetic Field 

In thi s work the IMF model originally intro- 
du ced bvlGiacalond and further developed 
by iGiacalone fc Jokipiil (|20oJ ). is used. Large- 
scale magnetic field fluctuations are generated by 
supergranular motion of the footpoints at the so- 
lar surface, with characteristic time T c and speed 

Vg. 



The resulting fluctuating IMF is given by: 



Dn 



Bprl V e {e,(j>,to) 
r 2 V w 
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where Vg(9, <f>, to) and V${0, 0, to) are velocities as- 
sociated with footpoint motion, V w is the solar 
wind speed, VLo is the solar rotation rate, Bo is 
the magnetic field magnitude at position r , r is 
radial distance from the Sun, 9 is colatitude, <p is 
longitude, t is time, and to represents the time at 
which a fluctuation at a given radius was produced 
on the solar surface. 

The velocity of the footpoint motions is mod- 
elled by a stream function ip(9, <j), to), satisfying: 



V e = - 



1 dtp 



sin 9 



Vs 



dtp 



(5) 
(6) 



The stream function is postulated to be: 

N n 



n—1 m— — n 



where: 

a™ = 6V g G™ (££ G n 



(7) 
(8) 



Gm 



, n 

1 <K 



10/3 



(-1/2) 



1 + «T C ) 5 / 3 



(9) 

A^c is the characteristic number of supergranu- 
lar cells, defined as N c = nro/V g T c , where r is 
taken to be one solar radius, r s . u>™ is the char- 
acteristic frequency of each mode and is taken as 
= 27m/I50X' c . The remaining parameters in 
the stream function are: random phase angles, f3™, 
drawn from a uniform distribution, and Y™{9,4>) 
the spherical harmonic function. The values of 
various parameters used in these simulations are 
presented in Table Q] 

The electric field is given by E = — ~V W x B/c 
since the plasma is quasineutral. 

In Section [3] the cartesian system used is de- 
fined as follows: the z axis is the rotational axis 
of the Sun, while the x and y coordinates lie in 
the heliospheric equatorial plane with the x axis 
corresponding to a longitude of <p = 0° . 

In the simulations presented here all particles 
are run through a realization of this IMF model. 
We set t = 0, giving to — —{r — ro)/V w , as the 
timescale of variation of the large scale turbu- 
lence is large compared to the particle propagation 
times considered. 

2.2. Particle Simulation 

The simulations presented in this paper use a 
full-orbit test-particle method, which involves nu- 
merically integrating the equations of motion of 
one particle at a time, repeating the process for a 
variety of initial conditions to assess the propaga- 



Parameter 


Value 


Bo (G) 
V w (cm s" 1 ) 
fio (rad s" 1 ) 

T c (day) 
V g (km s" 1 ) 
N 
ro (r s ) 


1.78 
5.0 x 10 7 
2.86533 x 1(T 6 
1 

2.0 
50 
1 



Table 1: Values of the parameters used in all the 
simulat ions. The values of T c , v g and N are those 
used in lPei et al.l ([20061 ) . 
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tion of a population. 

The code is fully relativistic and is a modi- 
fied version of one previously used to study par- 
ticle acceleration during magnetic reconnection 

The numerical tech- 
Bulirsch-Stoer 



(jDalla fc Browning! 12005 
nique us ed for integratio n is 



method (|Press et al.l TT9931 . The Bulirsch-Stoer 



routine is driven by a function which adaptively 
controls the stepsize of integration, so as to pro- 
duce results to a prescribed accuracy with min- 
imum computational effort. The solution's ac- 
curacy is determined by a user-specified tolerance 
value, which sets the fractional error allowed when 
integrating from t n to t n+ \. 

Extensive testing was carried out to find an ac- 
ceptable tolerance level for the solver; from these 
tests it was concluded that particle trajectories in 
the fluctuating IMF outlined in Section 12.11 con- 
verged when the tolerance was less than 10 , 
hence the simulations were carried out at a toler- 
ance of 10~ 12 . 

The presence of small-scale turbulence in the 
solar wind affects the trajectories of particles. 
Rather than providing a full description of the 
interaction of particles with turbulence in these 
simulations, its effect is described as a series of 
scattering events, each causing an instantaneous 
change in a particle's velocity. This 'ad-hoc' scat- 
tering ha s been used by m any others in the liter- 
ature (ee lPei et all (|2006l )V 

A mean free path, A, is fixed and the mean 
scattering time t scat is obtained by t sca t — A/u, 
where v is the particle's velocity in a frame of ref- 
erence which is stationary with respect to the solar 
wind. A series of scattering events are then numer- 
ically generated for each particle, with mean time 
tscat, in such a way as to form a Poisson distri- 
bution. A scattering event involves firstly trans- 
forming the particle's velocity vector into the solar 
wind frame, then the particle's velocity vector is 
reassigned by drawing a velocity vector at random 
from an isotropic distribution, changing both the 
particle's pitch- and phase-angles. The particle's 
new velocity is then converted back into the iner- 
tial frame of reference, and the integration of the 
particle's equations of motion continues. 



3. Results 

In all simulations presented in this paper a 
mono-energetic population of protons, with veloci- 
ties randomly distributed in a hemisphere directed 
away from the Sun, is injected instantaneously into 
the interplanetary medium. The location of injec- 
tion is a 3° x 3° area on a spherical surface of 
radius 21.5 r s centred on a colatitude of 60° and 
a longitude of 1.5°. 

3.1. Scatter- free Propagation 

Initial simulations were performed without ad- 
hoc scattering in order to investigate the cross- 
field transport of SEPs subjected only to large- 
scale fluctuations, and contrast their propagation 
in such a field with that in a Parker spiral field. 

Figure [T] shows the spatial distributions of 10 
MeV protons 30 and 60 minutes after injection. 
The left panel shows an x-y projection (along the 
heliospheric equatorial plane) and the right panel 
an x-z projection. Red dots are particles in the 
fluctuating IMF configuration and blue dots those 
in the Parker spiral (long dashed line). Field lines 
associated with the fluctuating IMF are indicated 
by short dashed lines. It can be seen that the lon- 
gitudinal and latitudinal spread of SEPs is much 
greater in the fluctuating IMF model than in the 
Parker spiral. 

The smallest scale of the fluctations is 27rr s /50 = 
87 Mm, this is large compared to the gyroradii of 
SEPs propagating close to 1 AU, hence the cross- 
field transport observed in the simulation results 
is simply due to the field line meandering and 
drifts. The particles propagating in the fluctuat- 
ing IMF model focus at the same rate as those 
in the Parker spiral, as can be seen in Figure [2j 
showing the evolution of the mean pitch angle of 
the populations with time. 

To investigate any dependence of SEP propa- 
gation characteristics in the fluctuating IMF on 
particle energy, proton energies were varied from 
500 KeV to 1 GeV. Results showed that SEP pop- 
ulations followed the same path through a single 
realization of the fluctuating IMF regardless of en- 
ergy, and hence cross-field transport in this con- 
figuration is due solely to field line meandering. 
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Fig. 1. — An x-y (left) and x-z (right) projection of 10 MeV proton locations at £=30 mins and i=60 mins 
after injection for scatter- free propagation in fluctuating IMF (red dots) and Parker IMF (blue dots). The 
thick dashed line is a Parker spiral line from the centre of the injection region, and dotted lines are field lines 
of the fluctuating IMF with starting locations at the center and corners of the injection region. 



3.2. Propagation with Scattering 




1000 2000 

Time (Sec) 



Fig. 2. — A plot of average pitch-angle vs time for 
all particles in the Parker spiral (dashed line) and 
the fluctuating IMF (solid line) simulations. The 
proton energy is 10 MeV and the distance from 
the Sun at the final time is about 215 r s . 



The scatter-free simulation described in Sec- 
tion 13.11 showed that in such conditions cross- field 
transport in the fluctuating IMF is due solely to 
field-line meandering. In this Section the effect 
of ad-hoc pitch-angle scattering, incorporated into 
the model as outlined in Section [2~2l is studied. 

Figure [3] shows populations of 100-MeV protons 
propagating in the fluctuating IMF. Each column 
shows a particle population at three different times 
for a given propagation condition. The initial con- 
ditions are the same in the three cases. The cen- 
tral and right columns are affected by pitch angle 
scattering, with mean free paths of 2.0 AU and 
0.3 AU respectively, and the left column shows 
scatter-free propagation. In the scattered popula- 
tions there are a number of particles close to their 
injection point at late times, in addition to a num- 
ber of outlying particles that have strayed a large 
distance across the field compared to the scatter- 
free case. 

To quantify the cross-field transport within 
the simulations, the parallel and perpendicular 
displacements of the particles' location from the 
Parker spiral field line originating from their ini- 
tial positions were calculated, for runs of 10000 
protons of 50 MeV energy. The first step consists 
of deriving the so-called 'target point', i.e. the lo- 
cation on the Parker spiral field line starting at the 
particle's initial posi tion, located near est to the 
actual final position (jTautz et al. i |201lh (see Ap- 
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Fig. 3. — x-y projections of a population of 2000 100 MeV protons propagating in the fluctuating IMF 
at three different times. The left column is for scatter-free propagation, the central one for ad-hoc pitch- 
angle scattering with A = 2.0 AU, and the right one for scattering with A = 0.3 AU. Times indicated are 
dimensionless, but can be converted to hours using 1.0e + 7 = 1 hour. The grey lines are magnetic field lines 
with seed points on a grid overlaid on the injection region. For the A = 2.0 AU case, many particles have 
propagated to large radial distances and are not shown in the plot. 
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pendix POI and Fig. fTTj) . If r indicates the location 
of a particle at a given time and Yt is the location of 
the target point, the displacement of the particle 
with respect to its target location is As = r — r t . 
Next a local Parker spiral coordinate system, cen- 
tered in the target location, is introduced, with an 
axis pointing outwards along the spiral (the 'par- 
allel' direction), a second axis in the direction of 
eg>=— eg with eg the standard spherical coordi- 
nate system unit vector and a third axis e^/ com- 
pleting the orthogonal system (see Appendix lA.il) . 
The components of Asj_ in this system are As^' 
and Asgi . The distance travelled along the Parker 
spiral between the initial location and the target 
point is indicated as 

The distributions of displacements for the three 
propagation conditions of Figure [3] are examined 
first after a fixed time from injection (when the 
average radial distance from the Sun in the popu- 
lation will be different in the three cases) and sec- 
ond at the times when the average distance from 
the Sun is the same for the three populations and 
set to 1 AU. 

Figure [4] shows the distributions of the displace- 
ment's components for the scatter-free, A = 2 AU 
and A = 0.3 AU conditions, for 50 MeV protons, 
three hours after injection. The top panel, dis- 
playing the distribution of shows that the pop- 
ulation with A — 0.3 AU has traveled the shortest 
distance along the field, followed by the A = 2.0 
AU case, and that the unscattered case has trav- 
eled the furthest, as one would expect. This differ- 
ence in average field-parallel displacement affects 
the magnitude of perpendicular displacement, as 
the spread of the field lines in the perpendicular 
direction increases with distance along the field. 
This leads to the distributions in the middle and 
bottom panels of Figure HI where the scatter-free 
population has wide distributions since all these 
particles have been able to travel a greater dis- 
tance along the field. 

In an attempt to remove the effect of field line 
expansion on As± plots, Asg< and As^/ values 
were normalised by the scale lengths of a Parker 
spiral flux-tube at the corresponding distance from 
the Sun, as defined in Appendix IA.3I Figure [5] 
shows the normalized distributions. In the 9' com- 
ponent, the unscattered population still has the 
widest distribution, while in the cj>' component the 
widest distribution is in the A=2 AU case. It 




100 100 200 300 

As,, (rs) 



Fig. 4. — in (top), Asg> (middle) and Asp (bot- 
tom) distributions of 50 MeV protons at a time 
corresponding to around 3 hours. Unscattered 
population: light gray, A = 2.0 AU population: 
gray, and A = 0.3 AU population: black. 
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should be noted that the normalisation introduced 
only corrects for the expansion of the Parker spi- 
ral field lines: the superimposed large scale turbu- 
lence may be characterized by a different kind of 
expansion which is not accounted for here. 

Figure [6] compares As distributions at the time 
when each population has an average radial dis- 
tance from the Sun of 1 AU. For the scatter-free 
case, A = 2 AU and the A = 0.3 AU cases this is 
29, 33 and 91 minutes after injection respectively. 
Particles are grouped in this way in an attempt to 
find the most appropriate time to compare popula- 
tions with different scattering conditions. Popula- 
tions with little scattering will reach greater radial 
distances in a given time than populations with 
strong scattering, and at larger radial distances 
the field line expansion due to magnetic field fluc- 
tuations will be larger, artificially increasing the 
the Asj_ value assigned to a particle 

Grouping particles this way results in the A = 
0.3 AU population having the broadest distribu- 
tion in both perpendicular directions, however 
that population has travelled the furthest along 
the field (see In distribution in the top panel of Fig- 
ure [5]) and this produces the pronounced 'wings' 
in the middle and bottom panel of Figure [6] This 
effect is similar to that seen in Figure [4] and is due 
to field line expansion. 

Figure [7] shows the As± components normal- 
ized by the Parker spiral flux-tube scale lengths. 
The A = 0.3 AU population has the broadest 
distribution in both the </>' and 6' populations. 
The distributions for the scattered populations are 
broader than for the scatter-free one. 

An alternative way to remove the effect of field 
line expansion on As± values involves normaliz- 
ing Asj_ by Z|| for each particle in the popula- 
tion at a fixed time 3 hours after injection, giv- 
ing the cross-field displacement per unit distance 
along the field. The results of this normalization 
are presented in Figure [5J which shows that the 
population with A = 2.0 AU has the broadest dis- 
tribution of As^i flu and has many outliers in the 
plot of As0//Z||. This result ties in with Figure [3j 
where the A = 2.0 AU population qualitatively ap- 
pears to have outliers that have strayed the largest 
distance across the field. 

Particle transport can be characterized by dif- 
fusion coefficients, under the assumption that it is 



Scatter-free 

A = 2.0 AU 

A = 0.3 AU 




Fig. 5. — Ase> (top) and As^' (bottom) distri- 
butions of 50 MeV protons, normalized by Parker 
spiral flux tube widths 1$ and 1$ respectively, as 
defined in Section [A.31 at a time corresponding to 
around 3 hours. Propagation conditions are indi- 
cated as in Figure |H 
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Fig. 6. — in (top), Ase> (middle) and As^ (bot- 
tom) distributions of 50 MeV protons at the time 
when the average radius of each population is 1 
AU. Propagation conditions are indicated as in 
Figure |U 
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Fig. 7. — Ase> (top) and As^ (bottom) distri- 
butions for 50 MeV protons at the time when the 
average radius of each population is 1 AU, normal- 
ized by the Parker spiral flux tube widths defined 
in Section IA.3I Propagation conditions are indi- 
cated as in Figure |H 



9 





Scatter-free — 




A = 2.0 AU - 




A = 0.3 AU : 




1 




Fig. 8. — Ase'/Z|| (top) and As<p>/1\\ (bottom) 
distributions of 50 MeV protons at a time corre- 
sponding to around 3 hours. Propagation condi- 
tions are indicated as in Figure HJ 



of diffusive character. The values of diffusion co- 
efficients K0/ and Kg/ are calculated from the 
equation: 



((As||,^) 2 ) 
It 



(10) 



In the expanding magnetic field the propagation 
parallel to the field is a combination of diffusion 
and focusing-driven streaming. This can be repre- 
sented as advection of a diffusively-spreading par- 
ticle p opulation, in the strong scattering case ([Earl 
19761 ). In line with this concept, the parallel dis- 
placement As 1 1 is defined relative to the mean po- 
sition of the particle distribution as: 



As 



II - l \\ 



Cll) 



(11) 



where the average is taken over the particle pop- 
ulation. 

This definition of diffusion coefficient includes 
the contribution of displacements in the particle's 
trajectory due to large-scale fluctuations in the 
magnetic field, and as such may not be suitable 
to directly compare with the diffusion coefficient 
used in the Parker transport equation. 

Figure [9] shows the diffusion coefficients calcu- 
lated from Eq.QUlfor the three propagation condi- 
tions as a function of time. For diffusive behaviour 
it is expected that the diffusion coefficient would 
reach a constant value. The top panel shows km 
for the three populations, although the km values 
of the unscattered population are included only for 
reference, as this case is not diffusive since the par- 
ticles are rapidly focussed, and travel as a beam. 
The increase in K|| with time for the scatter-free 
case is due to field line meandering. 

The slopes of plots shown in Figure [9] can be 
characterized by the relative differences between 
values at the midpoint and final time point in 
the simulations, as given in Table [5J In the A = 
2.0 AU case Ak||/k||=0.28, whereas for A = 0.3 
AU Akm/k|i=— 0.16, so neither curve is constant. 
For the perpendicular components of the diffu- 
sion coefficient an approximately constant value 
is reached only in the A = 0.3 AU case. 

The ratios of various diffusion coefficients from 
each population are presented in Table [3l It can be 
seen that in both scattering conditions the value 
of k$i/k\\ is about 0.01, while k^'/km is approxi- 
mately 0.08 for A = 0.3 AU and 0.04 for A = 2.0 
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AU. The ratio of the two perpendicular compo- 
nents of the diffusion coefficient is not one, as 
Ke'/ K 0'^O-2- This difference may be due to the 
fact that elec tric field drift aids diff usion in the 
direction ( Burns fc Hafpern 1968f) . or may be 
a property of the magnetic field meandering. 

The increased value of k^i compared with Kg' 
found i n these simulations d iffers from the conclu- 
sion of I Jokipii et al. I (Il995l) . based on analysis of 
magnetic field fluctuations in Ulysses data. How- 
ever this difference may be explained by the fact 
that the Ulysses measurements were made at a 
much larger radial distance from the Sun. 

The diffusion coefficients were also calculated 
by means of values of Asgi and As^ normalised by 
the width of a Parker spiral flux tube, as outlined 
in Section IA.3| the plots of the resulting diffusion 
coefficients are shown in Figure QUI The observed 
decrease in the diffusion coefficients with time can 
be understood by noting that in a simple radially 
expanding field without scattering, field line ex- 
pansion would generate a dependence like t^ 1 for 
the normalised coefficients. 

The simulations were also run with two new 
sets of random phase-angles (/3™) in the stream 
function defined in Equation [3 and it was found 
that this does not produce a significant qualita- 
tive change in the results. Plots corresponding to 
figures H] - [TU] for the two additional realizations 
display the same trends as the figures shown in 
this paper. Actual values of the diffusion coeffi- 
cients are broadly consistent. However in one of 
the additional realizations the ratio Kq/k'^ at the 
final time was equal to 1.6 (while in the other addi- 
tional realization this value was 0.1). This appears 
to be related to a large deviation of the magnetic 
field lines, from the Parker spiral, in the 9' direc- 
tion. 

If the value of N is increased from N — 50 to 





A = 2.0 AU 


A = 0.3 AU 


A«||/k| 


0.284 


-0.163 


Ak^/k^ 


-0.506 


-0.080 


Ako'/ko' 


0.454 


-0.083 



Table 2: The difference between the midpoint and 
endpoint values of the diffusion coefficients, nor- 
malised by the final value, for different propaga- 
tion conditions. 





A = 2.0 AU 


A = 0.3 AU 


ko>/k\\ 


0.011 


0.0139 




0.043 


0.083 




0.259 


0.169 



Table 3: Diffusion coefficient ratios in different 
propagation conditions. The values of each dif- 
fusion coefficient at the final time were used to 
calculate these ratios. 




2.0x10 



Scotter-free - 
X = 2.0 AU - 
K = 0.3 AU - 



3 6.0X10 3 8.0X10 3 l.OxlO 4 1.2x10* 1 A 

Time (s) 




2.0x10* 4.0x10' 6.0X10 3 8.0xiO J 1.0x10* 1.2xi0 4 1.4x10' 

Time (s) 



Fig. 9. — K|| (top), ngi (middle), and (bottom) 
running diffusion coefficients for 50 MeV protons 
under different propagation conditions. Propaga- 
tion conditions arc indicated as in Figure 01 
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Fig. 10. — Kg' (top) and k^i (bottom) running dif- 
fusion coefficients for 50 MeV protons, normalized 
by the width of a Parker spiral flux tube, outlined 
in Appendix IA.31 Propagation conditions are in- 
dicated as in Figured] 



N = 75 it is found that the particles' trajectories 
are very similar to those in the N = 50 case, with 
some small differences due to the fact that there 
is less power in the low N modes. 

4. Discussion and Conclusions 

The effect of large-scale IMF fluctuations on the 
propagation of SEPs is investigated by means of a 
full-orbit test-particle method. In the presence of 
large-scale turbulence, energetic particle popula- 
tions experience more extensive cross-field trans- 
port than would be the case in a Parker spiral con- 
figura tion, due to field line meandering (jGiacalond 
1999(1 . 

In this paper, three propagation conditions 
were analysed and the cross-field transport in 
large scale turbulence characterized. In scatter- 
free propagation SEPs simply follow field lines and 
the perpendicular transport observed is the result 
of magnetic field characteristics. By comparing 
populations at the same average radial distance 
from the Sun, it was found that when pitch-angle 
scattering is present, the perpendicular transport 
is enhanced (see also Figure [3]). Qualitatively, 
the introduction of scattering produces outlying 
particles that travel much further across the mag- 
netic field than unscattered particles. The total 
number of outliers is small in these simulations, 
due to computing time constraints on the total 
number of particles considered, however they are 
representative of a population that can be found 
at large separation from the injection region. Ob- 
servations show that the SEP fluxes measured at 
large angular separation from the source can be 
orders of magnitude smaller than those detected 
by a well connected spacecraft. 

The introduction of scattering makes it possible 
for particles to jump onto nearby field lines, by a 
distance of the order of the Larmor radius, and, 
due to field line meandering, a superposition of 
a large number of these events can result in large 
cross-field displacement from the original field line. 

The Parker spiral geometry poses some chal- 
lenges to the interpretation of the simulations. In 
an attempt to resolve these challenges a local co- 
ordinate system was introduced to yield parallel 
and perpendicular displacements from the parti- 
cle's original Parker spiral field line. The coordi- 
nate transformation to this Parker spiral system 
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(as detailed in Appendix IA.1|) is generic and ap- 
plicable to other heliospheric problems. For differ- 
ent scattering conditions, particles have travelled 
different distances along the field, influencing the 
calculated cross-field transport due to expansion 
of the magnetic field. It is important to take this 
effect into account in quantifying cross-field diffu- 
sion. 

The method used here has been firstly, to anal- 
yse Asj_ distributions for different scattering con- 
ditions after a fixed time from injection. It was 
found that the population which undergoes least 
scattering has the broadest distribution of As± 
values: this is because the least scattered particles 
can reach locations far away from their injection 
point where expansion and field line meandering 
cause large As± values. 

Secondly, displacements were analysed at differ- 
ent times for the three regimes, with times chosen 
so that the average radial distance from the Sun 
of the population is 1 AU. It was found that the 
A=0.3 AU population has the broadest As± dis- 
tribution, however this may be due to their having 
a significant number of particles at large distances 
from the Sun. 

A normalisation of perpendicular displacements 
was introduced in an effort to minimise the effect 
of field line expansion. The normalisation only 
takes into account the Parker spiral expansion, 
thus the expansion of the large-scale turbulence 
may not be completely compensated. If the dis- 
tributions of Asj^/As|| are plotted, it is seen that 
the population with A = 2.0 AU has the most 
cross-field transport per unit parallel transport, 
while the unscattered case has the least amount 
of cross-field transport per unit parallel transport. 
This appears to fit well with the qualitative con- 
clusion reached from Figure |3] that more distant 
outliers are present when A=2 AU compared to 
A=0.3 AU. While it is clearly established that per- 
pendicular transport is enhanced in the scattering 
cases compared to the scatter-free one, the analy- 
sis is not conclusive as to which of the two scat- 
tering regimes results in the most efficient perpen- 
dicular diffusion, due to the challenge presented 
by field line expansion. 

The ratio of perpendicular to parallel diffusion 
coefficients varies from 0.01 to 0.08 for the condi- 
tions analysed. Diffusion across the field in the lat- 
itudinal direction appears slower than in the lon- 



gitudinal direction. This is consistent with SEP 
observations at hi gh heliolatitude from Ulysses 
(jDalla et alJl2003bl ). Transport in longitude could 
be aided by the electric field drift which is directed 
along e$. 
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A. Appendix Material 



A.l. Transformation to a local Parker Spiral coordinate system 



A new coordinate system with origin at a point (xt,yt,zt) in space and one axis coinciding with the 
Parker spiral magnetic field direction at the point is introduced. The e; axis points outwards along the local 
direction of the Parker spiral magnetic field, a second axis is in the direction of eg/—— eg with eg the standard 
spherical coordinate system unit vector and a third axis completes the orthogonal right-handed system. 
The latter axis does not coincide with the standard spherical coordinate system unit vector e^. 

The angle between the radial direction and the direction of the magnetic field at the point is indicated by 
ft and is given by: 



ft = tan 



B r 



= tan 



r flo sin 9 

Vw 



(Al) 



where the terms B r and are the radial and longitudinal components of the Parker spiral magnetic field, 
and can be found by setting Vg = = in Eqs. (QJ and ([3]). 

If a vector has components (v x ,v y ,v z ) in the cartesian coordinate system introduced in Section 1, its 



components (vi,v^i ,vgi) in the new Parker spiral coordinate system are given by: 



vi 
vg> 



sin 9 cos 4> cos ft + sin ft sin <f> sin 9 sin <j> cos ft — sin ft cos cos ft cos 9 
— cos ft sin 4> + sin 9 sin ft cos (f> cos ft cos <j> + sin 9 sin ft sin <j> cos 9 sin ft 
— cos 9 cos <b — cos 9 sin 6 sin 9 



(A2) 



A. 2. Determining the target point 

The target point is the location on the Parker spiral field li ne starting at the particle's initial position 



with the shortest distance to the particle's actual final position (jTautz et al.ll201lh (see Figure [TT]) 



The targ e t poin t is determined by means of an iterative procedure that differs from the one used by 



Tautz et al. ( 201ll ). As an initial guess for the distance from the Sun of the target location, rt, the value 



of the actual radial distance for the particle is used. For a given r t , the target longitude can then be found 
as 4>{rt) = 4>o — (^o( r * — r o)/Vw) given the particle's initial colatitude 9q, longitude 4>o, and initial radius 
ro. The target colatitude is 9q. If the cartesian components of this initial target location are indicated as 
x t , y t and z t , the cartesian components of the As vector with respect to this target can be obtained. Using 
the transform defined in Eq. (|A.1[) . this can be converted into its components in the local Parker spiral 
coordinate system defined in Appendix lA.il 

The optimal target location will have a zero component along the e; axis and can therefore be found by 
solving the equation: 

(sin 9 cos 4> cos ft + sin ft sin 4>) [x e —x t ] + (sin 9 sin <f> cos ft — sin ft cos [y e — yt] + (cos ft cos 9) [z e — z t ] = (A3) 
where (x e ,y e , z e ) is the particle's actual location. The equation is solved by means of Brent's method 



(Press aifflH) This ensures that the particle's displacement with respect to the optimal target location 



is perpendicular to the Parker spiral. 

Care must be taken when specifying the bracketing region to ensure that only one root of Equation IA3I 
lies within the region specified, and that this root corresponds to the closest target point to the particle end 
point. 

A. 3. Normalization by scale length of Parker spiral flux tube 

The cross-section of a Parker spiral flux tube increases with radial distance from the Sun. In this Section 
the scale lenghts characterising the cross-section in the <fi an d directions are determined, indicated as 
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and Lg respectively. These quantities are used to normalise perpendicular particle displacements so as to 
remove the effect of Parker spiral field line expansion. 

The borders of the particle injection region, located at a distance n from the Sun and with angular 
extent Atp and AO in longitude and latitude respectively, define a Parker spiral flux tube with cross-section 
approximated by a\—r\ A(f)Ad. If the cross-section of the same flux tube at a distance r from the Sun is 
indicated as a, the conservation of magnetic flux requires that o=o\B\j 'B, where Bi is the magnetic field 
magnitude at r\ and B its value at r. 

In the Parker spiral: 




(A4) 



where a — v sw /(Q sm0). Hence: 




(A5) 



giving the following expressions for the cross-sectional scale lengths: 




(A6) 



Lg = rA9 



(A7) 



L$ and Lg are calculated using n=21 



.5 r 8 and A^=A6»=7r/60. 
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Fig. 11. — An example particle trajectory (dotted) and the path of its original Parker spiral field line (solid) 
are shown, in the x-y plane (left) and the x-z plane (right). The particle's end point and target point are 
indicated by black dots. The target point was acquired by solving Equation IA31 hence the perpendicular 
displacement of the particle from its original Parker line is simply the distance between its end point, and 
its target point. 
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